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Abstract 

Precise asymptotics known for the Green's function of the Laplace 
operator have found their analogs for periodic elliptic operators of 
the second order at and below the bottom of the spectrum. Due to 
the band-gap structure of the spectra of such operators, the question 
arises whether similar results can be obtained near or at the edges of 
spectral gaps. As the result of this work shows, this is possible at a 
spectral edge when the dimension d > 3. 



1 Introduction 

Let 

d 

L{x, D) = Y, (D, + b,{x)) a,i{x) {Di + bi{x)) + c{x), (1) 

be a second-order differential operator in Mf^ with smooth real valued coeffi- 
cients ttji, bj, c. Here D = (Di, . . . , Dd) and Dj = —id/dxj. 
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The operator L is assumed to be elliptic, i.e., the matrix aj£ is symmetric 
and 

d 

aji^j^e >e\^\' (2) 

for some 6 > and any ^ = {^i, . . . ,^d) G 

We assume that all coefficients are Z'^-periodic where Z'^ is the integer 
lattice in Mf^, i.e., aji{x + n) = aji{x) for all x G W^,n G Z'^ and similarly for 
bj and c. 

If the operator is defined as in ([T]) on the subspace C^(W^) C L^(R'^) 
of compactly supported smooth functions, it is symmetric and essentially 
self-adjoint [33]. The corresponding self-adjoint operator L has the domain 

Most of the conditions above (second order, smoothness of the coefficients, 
lattice of periods being equal to Z"^, and sometimes even self-adjointness) can 
be relaxed, and we postpone a further discussion until Section HI 

In order to formulate the problem addressed in this text and state the 
results, we introduce some notation and notions from the spectral theory of 
periodic elliptic operators (see, e.g. [T ^I^TIIS^ ). 

The spectrum of the defined above (see ([T])) operator L in L^(]R'^) has a 
band-gap structure, i.e., it is the union of a sequence of closed bounded 
intervals (bands) [aj, C M, j = 1, 2, . . . such that aj, bj — )■ oo: 

a(L) = U[a„/3,]. (3) 
j 

The bands can (and do, when d > 1 [2S1EDIE1]) overlap, but they may leave 
open intervals in between, called spectral gaps. Existence of such gaps, in 
particular, is crucial for the properties of semi-conductors [3] and photonic 
crystals [16l[22]. 

Our interest here is in obtaining asymptotics of the Green's function for 
the operator L near and at the boundaries of its spectrum. These asymp- 
totics are relevant for many problems in theory of random walks, Martin 
boundaries, Anderson localization, and others (e.g., [2],|11[5],[271[28] ) . This has 
been done for the case of the points at and below the bottom of the spectrum 
of L [U [271 [28] (see also [36] for description of the results of [^ in discrete 
case) and at the internal edges in one-dimensional case [28]. The resulting 
asymptotics are similar to the ones for the Laplace operator —A, except the 
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appearance of oscillatory terms. While A = is the only boundary point 
of the spectrum of the Laplace operator, spectral gaps may appear in the 
periodic case, and thus the question arises of whether one can obtain precise 
information about asymptotics of the Green's function inside these gaps, or 
at their edges. The Combes-Thomas estimates (see [5] for their contempo- 
rary version) provide the decay estimates for the resolvent in terms of the 
distance to the spectrum. They, however, do not reflect the exact asymptotic 
behavior of the Green's function. It is the goal of this and the next paper to 
obtain asymptotics or (at least) estimates at and near the spectral gap edges. 
In this text, only the case of a spectral edge is addressed in dimension d >3, 
while the asymptotics of the Green's function inside a spectral gap (near its 
edge) will be considered elsewhere. The restriction d > 3 in the spectral edge 
case is natural, since at the bottom of the spectrum, the operator is critical 
in dimensions d <2 (see [26ll27ll3T] ). 

Let = [0, 1]^ C M'^ be the unit cube. W is a fundamental domain 
of R*^ with respect to the lattice F := Z^, i.e., W and its integer translates 
cover and overlap only on the boundaries. 

The dual (or reciprocal) lattic^llof F is F* := 27rZ'^ and its fundamental 
domain is B = [— 7r,7r]'^. We call B the Brillouin zone (although in the 
solid state physics usually a somewhat different fundamental domain of the 
reciprocal lattice bears this name [3]). 

We also use the notation T = M'^/Z'^ and T* = R'^/27rZ'^ for the d- 
dimensional tori corresponding to the lattices Z"^ and 27rZ'^, respectively. 

Definition 1.1. For any k G C"', the subspace = H^iW) C H'{W) 
consists of the restrictions to W of functions f G Hf^^{M.'^) that satisfy the 
Floquet-Bloch condition that for any 7 G Z'^ 

/(x + 7) = e^'-V(a;) a.e. (4) 
Here denotes the standard Sobolev space of order s. 

Note that = L^{W). 

We adopt the name quasimoments(§ for vectors /c in (jlj). 

The periodicity of the coefficients of the differential operator L{x, D) 
preserves condition (jl]). It thus defines an operator L{k) in L'^{W) with 

^For a lattice F, the elements k of the dual lattice are characterized by the property 
that K • 7 G 27rZ for any 7 e F. 

^The name comes from the solid state physics [31132]. 
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domain Hl{W). We could alternatively have defined L{k) as the operator 
D + k) in L^(T) with domain H^{T). In the first model, L{k) is realized 
as the /c-independent differential expression L{x, D) acting on functions in 
W with boundary conditions depending on k (which can be identified with 
sections of a linear bundle over the torus T), while in the second one, the 
/c-dependent differential expression L{x, D + k) acts on the /^-independent 
domain of periodic functions on W. 

Notice that condition (j4]) does not change when k is modified by adding 
an element of l-n'L^ . This, when dealing with real values of quasimomentum 
fc, considering vectors k & B is sufficient. 

The hypotheses on L force the spectrum a{L{k)) of the operators L{k), 
k G W^, to be discrete. The following classical result (see and the 

references therein) describes the spectrum of the operator L in L/^iW^): 

Theorem 1. The following equality holds between the spectra of L and L{k): 

a{L) = U aim). (5) 

In other words, the spectrum of L is the range of the multiple valued 
function 

k ^ X{k) := a{L{k)),k e B. (6) 

The graph of this function in i? x M is called the dispersion relation or 
dispersion curv^ for the operator L: 

{{k,X):\ea{L{k))}. (7) 

Taking into account that the operator L{k) is bounded from below, one can 
label its eigenvalues in non-decreasing order: 

Xi{k) < X2{k) < . . . . (8) 

One thus gets the single-valued, continuous, and piecewise analytic band 
functions Xj{k). The ranges of these functions constitute exactly the bands 
of the spectrum of L shown in . 

In this text, we study the Green's function at an edge of the spectrum. 
Without loss of generality, by shifting the operator by a constant and chang- 
ing its sign if necessary, we can assume that the edge of interest is and the 

■^Another name of this graph is Bloch variety of the operator L, e.g. |21j . 
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adjacent spectral band is [0, a] for some a > and there is no spectrum for 
small values of A below zero. Thus, there is a spectral gap below zero. This 
means in particular that zero is the minimal value of at least one of the band 
functions Xj{k). 

If zero is the bottom of the entire spectrum, then in many cases (e.g., for 
periodic non- magnetic Schrodinger operators), it is known that the minimum 
is attained by a single band Xj{k) function at a single value of the quasimo- 
mentum ko (modulo 2-7rZ^-periodicity) and is non-degenerate, i.e., the Hes- 
sian of \{k) is non-degenerate at ko p!0] - [T2| [T7 tl3lfi Although this is not 
necessarily always true for other spectral edges, it is expected and commonly 
assumed in the mathematics and physics literature that it is "generically" 
true, i.e., it can be achieved by small perturbation of coefficients of the oper- 
ator. We thus impose the following assumption, under which we will be able 
to establish our main result: 

Assumption A 
There exist k^ E B and a hand function Xj{k) such that: 

Al \j{ko) = 0, 

A2 \\,ik)\ > 6 > fort ^ J , 

A3 ko is the only (up to 2711? -periodicity) minimum of Xj, 

A4 Xj{k) is a Morse function near ko, that is, its Hessian H at ko is a 
positively definite quadratic form. 

Thus, in particular, 

Hk) = \(k- kofH{k - ko) + 0{\k - ko\^) (9) 

is the Taylor expansion of Xj at ko. 

The parts Al and A2 of this assumption are known to be generically true 
for Schrodinger operators [18]. It is commonly believed, although not proven, 
that conditions A3 and A4 also are generically satisfied. Our result can be 
easily reformulated for a somewhat weaker condition than A3. Namely, one 
can assume 

A3' There are finitely many (up to 21x7? -periodicity) minima of Xj. 
^In the non-magnetic Schrodinger case, ko = \17\ . 
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Validity of the non-degeneracy condition A4 is often assumed, for instance 
in order to define the effective masses [3] or when studying emergence of 
impurity states under locahzed perturbations of the periodic medium [HHSj. 

In what follows, the branch Xj{k) plays a special role, and from now on 
we use the shorthand notation 

A(A;) := Xj{k). (10) 

Let us also introduce the notation (f{ko,x) for a Z'^-periodic function of 
X such that ip{x) = e'^^'^ipiko.x) is the (unique up to a constant factor) 
eigenf unction of L{kQ) with the eigenvalue A(fco)- 

L{ko)^ = X{ko)^. (11) 

We are ready now to state the main result of this text: 

Theorem 2. Let d > 3, the operator L satisfy the assumption A, and = 
{L + e)""*^ for a small e > denote the resolvent of L near the spectral edge 
A = (which exists, due to Assumption A). Then: 



1. For any (f),ij e L^^rnpi 



R-e(f) ■ 1p R4> ■ 1p 



for an operator R : Ll^^p(R'^) i-)- Lf^^{R'^). 

2. The Schwartz kernel G(x, y) of R, which we will simply call the Green's 
function of L, has the following asymptotics when \x — y\ ^ oo: 

G{x,y) 

= vr-^/^r(^)e-(-^)-^° ^(fco,x)Rfc^ (l + 0(\x- v\-')) 
2y/d^\H-y\x-y)Y-^ 11^(^1112^ ^ ^ 

(12) 

+ r(x,?/), 

where r(x, y) = 0{\x — y\~'''^^) ■ 
Here H is the Hessian matrix from ([9]) and the standard notation Ll^p{W^) 



and L^Q^(M'^) are used for the spaces of square integrable functions with com- 
pact support and locally square integrable functions on correspondingly. 
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The idea of the proof is to show that only one branch of the dispersion re- 
lation (namely, the one appearing in Assumption A) governs the asymptotics 
with the crucial term being the Hessian H at the point k^. This appears like 
an analog of homogenization (which is an effect that occurs at the bottom of 
the spectrum only), and indeed it is a version of it, along with other effects 
described in [HlIlillSllSMS]. 

We prove Theorem [2] in Section [3] after introducing the tools of Floquet 
theory [2T1[32] in Section |2l Section |4] contains final remarks and is followed 
by an Acknowledgments Section. 

2 A Floquet reduction of the problem 

In this section we use Assumption A and reduce the problem of obtaining 
asymptotics for the Green's function G(x, y) to the problem of finding asymp- 
totics of a scalar integral expression. The scalar integral we obtain resembles 
the Green's function of the Laplacian. 

2.1 Floquet transform 

Let 7 G M'^ and a function /(x) be defined on M"'. We then denote by the 7- 
shifted version of /: f-y^x) = /(x+7). Assuming that f{x) decays sufficiently 
fast (e.g., / is compactly supported, although f & would suffice, since it 
guarantees the L^^^^-convergence in (fT4|) below), the following transform plays 
the role of Fourier transform when studying periodic problems (in fact, it is 
the Fourier transform with respect to the lattice T/' of periods): 

Definition 2.1. The Floquet transform (also sometimes called Gelfand 
transform ) 

/(x)^/(fc,x) (13) 

maps a function f{x) 0/ x G M'^ into a function f{k,x) of {k,x) E B x W , 
defined as follows: 

f{k) = f{k, x) = j2 M^K'"-' = E + ^)^"'''- (14) 

Notice that f{k,x) is 27rZ'^-periodic with respect to k (recall that F* : = 
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2t\U^ is the dual lattice to F := Thus, it can be naturally interpreteci^ 
as a function on the torus T* with respect to k. With respect to the x- 
variable, though, the function is not periodic, but rather satisfies the cyclic 
(or Floquet) condition 

/>,x + 7) = e^'V>'^) (15) 

for all 7 G Z'^. Thus, the function f{k,x) is completely determined by its 
values on B X W. We usually consider it as a function f{k,-) of k G B with 
values in a space of functions defined on W. 

The next result collects some well known properties of the Floquet trans- 
form (e.g., pi]). 

Lemma 2.2. 

1. The Floquet transform is an isometry (up to a scalar factor) between 
L2(M'^) andjL^iW)dk = L^iB,L\W)). 

B 

2. Its inversion for any f G L^(R") is given by the following two equivalent 
formulas: 

f{x) = j f{k,x)d'k, xeW^ (16) 

and 

f(x) = j f{k, X - -i)e'^-^d'k, xeW + (17) 

3. Under the Floquet transform, the action of the operator L in L^(R'^) is 
transformed into 

e 

B 

i.e., 

rf{k) = mm 

for any f G H^{R'^). 

^Formally speaking, this interpretation requires to change the quasimomentum variable 
k to the Floquet multiplier variable z = {e^^^ , ...,e*'^''). We will abuse notation (e.g., in 
Lemma l2.2p . staying with the same variable k. 



8 



4- Let (j){k,x) be a function on M.f x such that it is 2ti'L'^ -periodic with 
respect to k, belongs to the space for each k (see Definition 
and is an infinitely differentiable function of k with values in H^. Then 
(j) is the Floquet transform of a function f G H^{M.'^). Moreover, for 
any compact set K C W'- the norm \\f\\H'>(K+'y) decays faster than any 
power of |7|. In particular, if s > d/2, then 

|/(x)| <C^(l + |xr^/or any N. 

Here d'k denotes the normalized Haar measure on T*, or the correspond- 
ing Lebesgue measure on B. 

2.2 The Floquet reduction 

We start by considering the resolvent in the spectral gap: 

R.J = {L + e)~V for / ^ L^^mpi^'^) and a small positive e. 
Applying the Floquet transform, 

lCf{k) = {m + e)-'f{k). 

Given another compactly supported function ip G L^^^piM.'^), we begin our 
analysis with the sesquilinear form 

(R-J^^P) = j^^ [im + e)"V>),^(A;)) d'k. (18) 

In the following lemma, we prove the first claim of Theorem [2] and introduce 
an expression for the operator R = lime_i.o R~t that is the basis of our study 
of the asymptotics of the Green's function. 

Lemma 2.3. Let the Assumption A be satisfied. Then the following equality 
holds: 

lim(i?_J,^) = (l(A;)"V>),^(A;)) d'k. (19) 
Thus, the Green's function is the kernel of the operator R defined as follows: 

Rf{k) = L{kr'f{k). (20) 
The integral in the right hand side of / fTPj) is absolutely convergent for fjipE 

[2 

comp ' 
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Proof. The complex Bloch variety of the operator L, denoted by S, con- 
sists of all complex pairs {k, A) G C^"*"^ such that the operator L{k) — \ on the 
torus T is not invertible. This variety is an analytic subset of co-dimension 
1 (e.g., |21]). Moreover, there exist an entire scalar function h{k, A) on <C^^^ 
and an entire operator- valued function I{k, A), such that 

1) h vanishes only on (the nonsingular) S, has simple zeros at all points 
of S, and 

2) outside S, the inverse operator {L{k)-\)~^ = h{k, X)~^Iik, A) [211135]. 
In particular, for a small e > 0, our operator of interest {L{k) + e)~^ = 

h{k, —e)~^I{k, — e). Due to assumption A, for small e > and quasimomenta 
k with small imaginary parts, the function h{k, — e)~^ is equal, up to a non- 
vanishing analytic factor, to {X{k) -|- e)~^. Thus, for such values of k and e, 
one can write the bilinear form {R^^f, (p) as 

T* 

Here M{k, e) is an operator-function in L'^{T), analytic for small e and nearly 
real k. The function A(A;), due to Assumption A, is a non-negative function 
with a single simple non-degenerate zero on T. Also, due to the compactness 
of support of f{x) and ipi^) and (fT4l) . the functions f{k),ip{k) are analytic 
with respect to k. Thus, for e > 0, the integrand in ( 12T|) is nonsingular. Since 
d> 3, it is still integrable when e = 0. Hence, a straightforward application 
of the dominated convergence theorem finishes the proof. □ 

Remark 2.4. Lemma E3 could be proven without using the results of 121^ 35 ^ 
concerning the representation of the resolvent {L{k) — X)'^ = h{k, X)~^I{k, A) 
(which requires some regularized infinite determinant techniques). One could 
instead achieve the same using eigenfunction expansions that are employed 
below. 

2.3 Singling out the principal term in R 

As we have just observed, the Green's function in question is the kernel of 
the operator R that acts according to fEUl) : for a compactly supported L^- 
function /, the Floquet transform of Rf is given by Rf{k) = L{k)~^ f{k). In 
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other words, 

Rf(x) = j L{ky^f{k,x)d'k for X e (22) 

or 

i?/(x + 7) = j e'^-^L{k)-^f{k,x)d'k for X G iy,7 G Z'^. (23) 

T* 

In this section, we will single out the part of the integral in (123|) that is 
responsible for the principal term of the Green's function asymptotics. The 
first thing to notice is that in order to determine the Schwartz kernel of R, it 
is sufficient to apply the operator to smooth functions. Thus, we may assume 
that / G C^iW^). 

Let us localize the integral around the minimum point k^ of \{k). In 
order to do so, we first recall [2HI35] that in a neighborhood of /cq, there exist 
functions %lj{k,x) and ip{k,x) that are analytic in k and satisfy 

1. ip{k,x) = e^^'^ipik.x), 

2. ip{k, x) is Z'^-periodic in x, and 

3. L{k)^{k) = X{k)^{k). 

There exists an analytic spectral projector P{k) that projects L^(T) onto the 
span of ipi^k). 

Let z^, ?7 G C^{B) (here B = [— 7r,7r]'^) be cut-off functions so that u = 1 
near ko, ?7|supp!^ = 1, and ip{k) is well-defined on supp?]. We can decompose 

f = uf + {l-u)f. 

The operator L{k) is invertible (in x) whenever k stays away from k^. Thus, 

Ug{k,x) = L-\k){l - u{k))f{k,x) 

is well-defined and smooth with respect to {k, x) (recall that the operator 
L{k) is elliptic and depends analytically on k). The smoothness of Ug means 
(according to the statement 4 of Lemma 12.21) that Ug has rapid decay in x 
(we thus chose the subscript g in Ug to indicate that this is a "good" part of 
u). 

Next, we must solve 

L{k)u{k) = iy{k)f{k). 
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We will look for a solution u of the form u = ui + U2 where Ui = P{k)ui, 
U2 = Q{k)u2, and Q{k) = J — P{k). In this spirit, let /i = P{k)u{k)f , and 
/2 = Q{k)i'{k) f . Observe that 



h = P{k)rj{k)u{k)f = r^{k)P{k)v{k)f = r^{k)f_ 



since 77 is a function of k and not of x. With the two orthogonal projection op- 
erators P{k) and Q{k), the problem of solving L{k)u{k) = f{k) decomposes 
into 

f fl 



L{k)P{k) 








L{k)Q{k) 



In this block-matrix structure, the first column corresponds to a one- dimensional 
subspace (the range of P{k)), while the second one corresponds to the infi- 
nite dimensional range of Q{k). Since the operator L{k)Q{k), considered on 
the range of Q{k), has no spectrum near zero, the operator {L{k)Q{k))~^ is 
well-defined and so is U2- Additionally, U2 is smooth, so according to Lemma 
12.21 again. U2 decays rapidly as |x| — )> 00. 

We know that iti = P{k)ui and L{k)ip{k) = \{k)ip{k). Putting these 
facts together, we need iii to satisfy the equality 

L{k)P{k)u, = P{k)h = v{k)P{k)h. 
which can be rewritten consecutively in the following forms: 



Therefore, it is useful to define 

. /, X _ v{k){fi{k,-),'^{k,-))Tilj{k,x) 
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Applying the inverse Floquet transform, we have for n G Z'^, 

2.4 From ui to a "reduced" Green's function 

We now pass from ui{x + n) to computing, up to lower order terms at infinity, 
the Green's function. Below, we will define a function Gq{x, y) that provides 
the leading term of the asymptotics of the Green's function G{x,y). To 
do so, our plan is to isolate a scalar integral, whose asymptotic behavior is 
responsible for the asymptotics of the Green's function. 
We introduce Go{x,y) as follows: 

ui{x + n)= Go{x + n,y)fi{y)dy. 

Observe that 



I n(k) V f e^"-^ L (.,).-rk., Hk,y-lMk,x) 



We notice that 

ij{k, y-i) = e^'-^y-^^ifik, y) = e-'^-^^{k, y) 



so ip{k, 2/ — 7) = e^^'^^ipik, y) and 



u.{x + n) = I My) I -i- / e-S(fc) f^^;"^,tiw!! ^fe 1 dy. 



Thus, we have 
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Note that if ?/ e T and 7 G Z"^, then 



or more simply, 



3 Asymptotics of the Green's function 

We start with deducing asymptotics in the reduced case. 

3.1 Reduced Green's function Go 

Let 

_ ip{k,y)if{k,x) 

Then p is a smooth function in k near k^, x, and y. Let fio{k — k^) = ri[k), 
so /io is a cutoff function supported near 0. Then 

It will be shown in the next subsection that the integral in (!25|) provides the 
main term of the asymptotics for the full Green's function G{x,y). 
By Taylor expanding around ko, we can write 

p{k) = p{ko) + {k- ko) ■ Pi{k). 

The integral then becomes 

We will see that the main term comes from the first integral on the righthand 
side of ( l26l) . It is this integral that we address first. 
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Assumption A2 means that the matrix H{X, ko) from Assumption A4 is 
positive definite. Let A = i7(A, /cq)- Then 

\{k) = ^{k - kofA{k - ko) +gi{k- ko), (27) 

where gi = 0{\k — ko\^). 

Let A^/^ be the positive definite square root of A. Then, after the changes 
of variables z = k — ko and = A^^'^z, we get 

i(x-y) k Mf^ - ^o) 

i{x-yyko f i{A-^l\x-y)YA^I^z [Iq{A-^I'^ A^l'^ z) 

\{Ay^zYAyH + gM-^'''A^'^^) 



{det Ay/^J^, ^e + ^?(0 

Here g = 2gio A~^/^ = 0{\^\^) and = /^o o A"^/\ 

Let Xo = 74~^/^(x — y). We have now reduced the problem to finding the 
asymptotics of the following integral in M": 

where g{C,) = 0{\^\'^). The main goal is to eliminate the extra term g(^), 
after which we have a situation similar to the ordinary Laplacian. To do 
this, we prove the following proposition: 

Proposition 3.1. Let A, g, and fi be as in j^) . // |xo| > 2, then 
where 

In the following proof we use the notation for a constant that depends 
on z only. The value of may change from line to line. 
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Proof. If the support of /i is small enough, then since g{^) = 0{\^f), we have 
where ri{C,) G Cq(M'^). Next, by Taylor expanding g{^), we can write 

E fiS-'-«) (30) 

^^^^ n=l 3n<|a|<2n+2+d 

for some constants (7a,n and r G Cg(]R'^). 

We start by handling the lower order terms in ( l30|l first, namely the terms 
with and r{^). 

3.1.1 Terms of the form 

Lemma 3.2. Suppose that h{x) is such that its Fourier transform h{C,) is a 
function homogeneous of order —1. Then 

[ e"-«/i(0/i(0 < (31) 
for \xo\ > 2. 

Proof. Let us observe that 

h{x)jl{x — xo) dx = h{x — xo)jl{x) dx. 



Since /i G C(j"(M'^), its inverse Fourier transform /i is a Schwartz function and 
hence rapidly decreasing, i.e., there exist constants for £ > so that for 
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Next, for |xo| > 2, it follows that 
h{x — Xo)fi{x) dx 



< 



+ 



1 



I X~XQ I > 

x\>l 



\x — Xol'^ 
nd—l 



X\ 



dx 



dx 



< 



|x|<l 



d-1 



Cn dx 



Cd,t 



for some constant that depends on fi and d. This proves the claim of the 
lemma. □ 

We will now use an integration by parts argument to reduce to the case 
of degree —1 homogeneity. If a = (ai, . . . , a„), then 



a function that is homogeneous of degree |a| — 1 



Assume that xq 



[X, 



0' 



Xq) and l^ol ~ l^ol. Then for |a| > i—1 and q := \a\—i + l one has 



< 



c 



1||q|-^+1 



X, 



\a\-l+l 

E 

j=0 



de^\i\'^ Oil 



dx 



If j = \a\ then we can use (!3T|) with h{^) = {-mi) and compute 



that 



i\\a\-e+i 



1 



< c 



If j 7^ |a| — ^ + 1, then there exists c so that c ^ < |^| < c if ,^ G supp |^ 

since /i G and /i = 1 in a neighborhood of 0. In this case, we can integrate 
by parts as needed to obtain decay of 0(|xo|~'^'''^) (or of any higher power of 
|xo|). 

Combining our estimates, there exists a constant C = C{fj,,d,g) so that 



d+l 



=1 3n<\a\<2n+2+d 



< 



c 



\xo\'-'' 
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3.1.2 Estimate for the r(^) term 

Sinct 

that 



Since r(^) G Cq(]R'^), we may integrate by parts {d — 1) times and observe 



3.1.3 The main term estimate 



Cd 



\d-l • 



We will show that for the constant Cd from (129|) . the following standard 
asymptotics holds: 



(2vr) 



Cd 



d-2 



+ 0(|xo| 



-d+l\ 



Let 6' be a radial, smooth approximation of the identity with supp 6 C 5(0, |) 
and Jjgjj 6 dx = 1. Recall that the volume of the unit ball in M'^ is Ud = p^^^^^^^ 
and the Newtonian potential is 

1 1 



d{d — 2)uJd \x 



d~2 ■ 



Since \l/ is a harmonic function on B{xq, |), it follows from the mean- value 
property that 

'^*9{xo) = ^'(xo). 

Therefore, denoting by /3(x) the function whose Fourier transform is fJ^i^), 
we get 



d^ 



d^. 



The inverse Fourier transform applied to yields 



d^ = ^*9{xo) = ^{xo). 



Also, since 9{0) = /3(0) = 1, we can write 



ISP 



l<|a|<d+l 
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where fia is rapidly decreasing in |^| and /i' G Co(M'^). 

Now the argument provided earher in the proof leads to the desired error 
estimate 

/ ^^'"^(^ E rUO + ^^'{0)d^ = 0{\xo\'~'). (32) 

Finally, since /i^ depends on 9 and (3, so the 0-term in (l32l) is bounded 
independently of |xo|. □ 

In order to finish the proof of the asymptotics of the reduced Green's 
function Gq, it remains to bound the last term in (1261) : 



{27rr Xik) 

This integral can be estimated similarly to the technique used for controlling 
the error terms in Proposition 13. II Let us use a Taylor expansion in {k — ko) ■ 
Piik) around fco- In the notation from the proof of Proposition 13. we have 



d+l 



jk-ko) ■ pi{k) , . \^ ^ , T^fc\ 
—- fio{k-ko)=p{02^ 2^ + 

^ ' n=0 3n+l<|/3|<2n+2+d 

The remainder function R{C,) belongs to Cq (M'^), and as in the proof of Propo- 
sition 13.11 it follows that 

If we collect our estimates, we arrive to the following intermediate result: 

Theorem 3.3. Under the conditions of Theorem\^ the following asymptotics 
holds for the reduced Green's function: 



GJx v) - ! ei^-y>^^!^2]Mh^dk 



27r'^/2(detA)V2 {\A-^/^{x-y)\J ||^(Mlli.m ^ ' 



L2(T) 

+ 0(\x-y\-''^'] 
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3.2 The full Green's function asymptotics 

We now need to show that when descending from the true Green's function to 
the reduced one, we have not changed the principal term of the asymptotics. 

Let us recall that the restriction to the reduced Green's function Go oc- 
curred in Section 12.31 when we dropped the operator term 

e 

T := j T{k)d'k, 

T* 

where T{k) = i'{k){L{k)Q{k))~^ + (1 ~ ^^(^))(-^(^))~^ is a smooth operator 
valued function of k. We will show now that the Schwartz kernel of this 
operator in decays sufficiently fast to be included into the error term 
r(a;, y) in f[T^ . 

To do so, we will need the following auxiliary result: 

Lemma 3.4. For any compactly supported functions (l),6 & C^lW^), the 
operator norm in L'^{M.'^) of the operator (p^iTO^, where g-y{x) := g{x + 7), 
satisfies the estimate 

||0yT^^|| = 0(|7-7'|"^) for any N > when -'y'l ^ 00. (33) 

Proof. By using unity partitions and shifts, one can assume without loss 
of generahty that supp (/) C W and supp 9 C W. Then the Floquet direct 
integral decomposition of the operator (f)yT9^ is (up to an absolute value 1 
exponential factor) 

e 

0(x) j e'''<^-^'^T{k)9{x)d'k, 

T* 

where the functions (p and 9 are considered as the multiplication operators by 
them in L^(T*). Then the estimate ([33D is just the decay estimate for Fourier 
coefficients of an infinite differentiable vector (operator) valued function on 
the torus (in this case, (f)T{k)9). □ 

The operator norm estimate (l33l) is weaker than the pointwise estimate 
of the Schwartz kernel r of the remainder in Theorem [2l However, one can 
bootstrap it from (133|) to these pointwise estimates. This can be done in 
different ways. One would be to use Schauder estimates, while another is to 
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apply the similar reasoning to the Floquet expansion of the operator T in 
a sufficiently smooth Sobolev space H^{M.'^) instead of L^, which would give 
corresponding smooth norms estimates on the Schwartz kernel of (p-y/TO^. 
This finishes the proof of Theorem [2l 

4 Final remarks 

• In the case when the weaker condition A3' is satisfied, one needs to add 
the asymptotics coming from all non-degenerate isolated extrema. 

• It is easy to reformulate (by spectral shift and, if necessary, changing 
the sign of the operator) the asymptotics without assuming that Aq = 
and that this is the upper, rather than lower bank of the gap. 

• One might wonder whether one truly needs to have the additive error 
r in f lT2|) of Theorem |2l not only the multiplicative one. This seems 
to be indeed necessary, unless one is dealing with the bottom of the 
spectrum. The reason is that the eigenfunctions (j) in f lT^ will have 
zeros, and thus the additive error r cannot be pulled in as a part of the 
multiplicative one. 

• Due to absence of restrictions on the coefficients of the operator L, the 
result holds for an arbitrary lattice of periods, not necessarily Z'^. 

• The result carries over without any essential changes in the arguments 
to periodic second order elliptic operators on abelian coverings of com- 
pact manifolds (as in [231426] ). 

• The smoothness conditions on the coefficients of L are exceedingly se- 
vere. One can easily follow the proof to detect that a finite smoothness 
would suffice. We decided not to pursue this line of study in this article. 

• As we have already briefly mentioned, the asymptotics we get has the 
spirit of an unconventional analog of "homogenization," which, unlike 
the classical homogenization, occurs at the internal gap edges, rather 
than at the bottom of the spectrum. One can find other results that 
can be interpreted as incarnations of such homogenization in pi [T^[TH| 
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